Measuring dark energy spatial inhomogeneity with supernova data 
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The gravitational lensing distortion of distant sources by the large-scale distribution of matter in 
the Universe has been extensively studied. In contrast, very little is known about the effects due to 
the large-scale distribution of dark energy. We discuss the use of Type la supernovae as probes of the 
spatial inhomogeneity and anisotropy of dark energy. We show that a shallow, almost all-sky survey 
can limit rms dark energy fluctuations at the horizon scale down to a fractional energy density of 
~ 10""*. 



Introduction — Type la supernovae (SNe) are thought 
to be excellent standard candles, with a calibrated lu- 
minosity distance known to ~ 15% By using these 
SNe to measure relative distances out to cosmological 
rcdshifts it is possible to reconstruct the expansion 
history of the Universe [H, H, Q . Over the past decade 
SNe have helped establish that the expansion of the 
Universe is accelerating [1, H, H, 0, Hi- Current and 
proposed SN surveys aim to significantly expand and 
refine the database of SN measurements, and thereby 
improve our understanding of the dark energy phe- 
nomena thought responsible for the accelerating ex- 
pansion. 

The next generation of SN observations will at- 
tempt to answer a fundamental question about dark 
energy: whether it is static, indicative of a cosmolog- 
ical constant, or time-evolving, such as quintessence 
and a wide range of alternative theories 0]. In ad- 
dition to variations in time, it is also prudent to con- 
sider the possibility of spatial variations. Fluctuations 
of the luminosity as a function of direction on the 
sky, relative to the global average, may indicate novel 
physics irrespective of the dark energy equation-of- 
state. The aim of this paper is to explore the degree 
to which future SN data may place interesting limits 
on the homogeneity and isotropy of dark energy. 

Dynamical dark energy is necessarily inhomoge- 
neous; covariance implies that a time-varying field is 
equivalently spatially- varying. Consequently, fluctu- 
ations in dynamical dark energy develop in response 
to inhomogeneities in the surrounding radiation and 
matter fields [lOl. [Til. [T^ . Where there are fluctuations 
in the dark energy density, there should also be spa- 
tial variations in the local expansion rate. These vari- 
ations may be directly observed as spatial anisotropy 
of the luminosity distance-redshift relation. In ad- 
dition, modifications to gravity can masquerade as 
dark energy (e.g. [i3l); spatial fluctuations in the 
new gravitational fields can leave a similar imprint 
on the local expansion rate and luminosity distances, 
and there are ongoing efforts to use fluctuation spectra 
to test for general departures from General Relativity 
[l3.[Tl.[l6l.1l7l.[T8t. Ahhough we do not expect spatial 
fluctuations in slowly-evolving quintessence flelds, or 
similarly weak modifications of gravity, to have mea- 



surable impact on the luminosity distance-redshift re- 
lation, our lack of understanding of the physics under- 
lying the dark energy suggests that we not preclude 
observable spatial fluctuations. 

Dark energy has also been speculatively tied to the 
breakdown of two fundamental tenets of cosmology, 
homogeneity and isotropy. If the observed cosmic ac- 
celeration is due to gravitational effects in a strongly 
inhomogeneous universe [Tol . [20| then we may expect 
similar distortions of the expansion rate and lumi- 
nosity distance. Anisotropic dark energy 
may be connected to a breakdown of global, statis- 
tical isotr opy , hints of which have been seen in the 
CMB [23, [131 . It is conceivable that these features are 
also present in the expansion rate of the Universe, and 
will be reflected in residuals in the SN luminosity dis- 
tance [2^ . Furthermore, tests for angular variations in 
the luminosity distance can be used to constrain novel 
photon interactions, such as photon-to-axion conver- 
sion [23, [2^ , which can distort our view of the cosmic 
expansion history. 

While we do not yet have an understanding of the 
physical nature of dark energy, certain candidate mod- 
els predict spatial fluctuations in the energy density of 
dark energy, relative to the global average, with ampli- 
tude 10~^ on horizon-size length scales This 
scale suggests that the best probe of dark energy clus- 
tering may be the CMB, in combination with other 
measures of large scale structure. An all-sky CMB 
map can be used to study dark energy anisotropics 
based on differences in the distance to the last scat- 
tering surface in different directions. Unfortunately, 
the CMB integrates out to z 1100, and thus is in- 
sensitive to dark energy physics at relatively low red- 
shifts. CMB limits for dark energy horizon-scale rms 
fractional fluctuations are not expected to be better 
than 10"^ 

The clustering of dark energy in a supernova sam- 
ple manifests itself as correlated fluctuations of the 
observed luminosity residuals. Of course, a variety of 
conventional gravitational phenomena also distort the 
paths of photons from distant SNe [s^, [3l[ . The dom- 
inant effect at small angular scales (< 5 arcmin) and 
high redshift (z > 0.5) is the lensing magnification due 
to the intervening mass fluctuations [U, [H, [s^ . This 
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lensing effect results in correlated spatial anisotropics 
in the luminosity distance estimated from nearby SNe 
on the sky, and can be measured in terms of magni- 
tude residuals relative to the best-fit cosmological dis- 
tance relation at a given redshift [s^ . Previous stud- 
ies have shown that a deep SN survey can be used to 
measure magnification statistics of the Universe (36j . 

There are also correlated fluctuations due to Sachs- 
Wolfe and integrated Sachs- Wolfe- like effects, consist- 
ing of weighted integrals of time derivatives of the 
scalar gravitational potential along the line of sight, 
but these are expected to be too small to be relevant 
[sot - Finally, there are spatial correlations due to the 
bulk flows in which the SN host galaxies are embed- 
ded [13, m, [111 , although these are only important at 
lower redshifts {z < 0.2). Efforts to discover corre- 
lated luminosity fluctuations induced by dark energy 
must overcome these nuisance effects. 

C'alculational Method — We calculate the fluctuations 
in the apparent luminosity due to perturbations in the 
dark energy component [4^ . We also incorporate the 
effects due to inhomogeneities in the matter, which ac- 
count for conventional gravitational lensing. We lack 
a theoretically compelling model for the dark energy, 
and thus confine ourselves to a purely phenomenolog- 
ical description of dark energy inhomogeneities, moti- 
vating a model-independent measurement that can be 
used to guide future tests of specific models of clus- 
tering. 

We model the luminosity distance, dL{z,h), to a 
source at redshift z as a function of sky location n. In 
a general Robertson- Walker space-time, the luminos- 
ity distance along a given line-of-sight is given by 
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where c is the speed of light, Hq is the Hubble con- 
stant today, and Qm, said arc the matter, cur- 
vature, and dark energy density in units of the crit- 
ical density. We specialize to the spatially-flat case 
for this study. We note that for the line-of-sight in- 
tegral, the relevant densities are the actual densities 
in the given direction (including the effect of spatial 
variations), not the globally-averaged densities. The 
apparent luminosity of a source is given by L{z,h) ~ 
Lx {dLldL(z,i\)Y where L and d^ are the mean lumi- 
nosity and luminosity distance, respectively, defined 
by the underlying, unperturbed cosmology [3j|. Per- 
turbations are given by (5L(z,n) = [L{z,ii) — L]/L. 
We build a phenomenological model for the dark en- 
ergy inhomogeneities, characterizing the fluctuations 
as 5x{z, n) = [ria;(2;, n) — ^x{z)\/?lx{z), relative to the 
mean dark energy density, at each redshift. We 
can then write the three dimensional power spectrum 
for the dark energy fluctuations in Fourier space as 
(<5,(k)5,(k')) = (2^)25(k + k')P.(fc). By observing 
these perturbations as a function of angular scale and 



redshift, we can attempt to quantify the underlying 
inhomogeneity of the dark energy [40| . 

In the limit of small number statistics, as would be 
the case with current and next-generation large SN 
surveys, it may be more practical to measure the an- 
gular power spectrum of luminosity fluctuations. As- 
suming statistical isotropy, the angular power spec- 
trum is given by {{5 L)\rai^ l) i' m') = Cf^ 8u,5„„rL' , 
where {5L)im arc the multipole moments of the lu- 
minosity fluctuations. The expected angular power 
spectrum of (fractional) luminosity fluctuations, 5l, 
for a sample of SNe at redshift z, is given by: 
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Here c?a(z) is the comoving angular diameter distance, 
r is the comoving radial distance, and n{r) is the radial 
distribution of SNe. 

Equation 2 neglects the contributions of gravita- 
tional lensing due to matter inhomogeneities. In what 
follows we will consider lensing as a source of noise, 
such that any measurement of Px{k) above this noise 
represents a measurement of the dark energy inhomo- 
geneity power spectrum. Instead of a specific model 
[4l|, we characterize the power spectrum of dark en- 
ergy fluctuations with a phenomenological model in 
Fourier space: 
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where 5xH is the amplitude at horizon scales, Ux is the 
tilt to the spectrum, s is the exponent in the growth 
factor, and kc is the cut-off scale for fluctuations. For 
quintessence fields the cut-off wavenumber is set by 
the dynamics of the equation of state. Formally, Px is 
the power spectrum of an additional clustering com- 
ponent to potential fluctuations, which we evaluate in 
the conformal-Newtonian/longitudinal gau ge ( we ab- 
sorb a prefactor proportional to ^x to 5xh [43 )• The 
parameters {SxH, nx, s, kc) can be constrained with 
future SN data. Note that while we have written a 
power spectrum associated with dark energy fluctua- 
tions, the same power spectrum may also describe the 
impact of modiflcations to gravity, or exotic physics 
of the dark sector. 

Gravitational lensing due to the inhomogeneous 
matter distribution sets a lower limit on the strength 
of detectable dark energy fluctuations. The dark en- 
ergy signal must rise above this amplitude, and thus 
the la upper limit on SxH (given other parameters) 
of a null experiment is given by 
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where 
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Here is the anisotropy power spectrum of SN hi- 
minosities expected from cosmic magnification due to 
mass fluctuations [s^, Nqn is the surface density of 
SNe (number per steradian), ctl is the uncertainty in 
the luminosity measurement from each SN, /sky is the 
fraction of sky covered by the survey, and A£ is the 
binning width in multipole space. When considering 
the extent to which finite SN surveys can search for 
inhomogeneities, we take Ij^in = 180° /9 where 9 is the 
one-sided area of the survey. We assume an uncer- 
tainty in the (calibrated) intrinsic luminosity of SNe 
given by ctl = 0.1. 

Results — Fig. [T] summarizes our results, illustrating 
the sensitivity of future SN surveys to clustering of 
the dark energy. We show the angular clustering 
power spectrum of luminosity variations relative to 
the mean luminosity. The weak Icnsing power spec- 
trum (thin solid curve) dominates at arcminute an- 
gular scales and at high redshifts We also show 
examples of dark energy clustering, arising from mod- 
els with = 10"^ s = 2, kc = 10-1 h Mpc'^ for 
three values of Ux ■ These spectra are chosen such that 
fluctuations in the luminosity residuals become impor- 
tant at large angular scales, and thus dominate over 
the perturbations due to matter fluctuations (gravi- 
tational lensing). For comparison, the amplitude S'^j^ 
for horizon-scale perturbations, relative to the dark 
energy density, is of order ^ 10~^ for certain dark 
energy models [U, [2§| . 

For comparison, in Fig. [1] we also plot the shot-noise 
for two potential SN surveys (dashed lines): 3,000 SNe 
in 10 deg^, as can be achieved with deep repeated ob- 
servations from space, and 10^ SNe in 20,000 deg^, 
from a shallow, but wide, SN survey. The latter is 
equivalent to what can be achieved with the Large 
Synoptic Survey Telescope (LSST; [13]). We combine 
the effects of lensing and shot-noise, generating a total 
noise estimate in wide multipole bins, following equa- 
tion ([6]) . The stair-step curves show this binned noise 
for the two SN samples we consider, where we take 
A£ = U for the bin sizes. 

Using the procedure outlined in equations ^ and 
([6]), we have estimated the minimum rms amplitude of 
dark energy spatial fluctuations at horizon scales that 
can be detected with future SN samples. A survey of 
3,000 SNe in 10 deg^ (e.g., from a space-based JDEM 
[4J]) provides limited constraints due to the narrow 
sky coverage, finding limits of 5xH < 0.2 if n^. —3, 
and SxH < 10"* if Ux ~ 0. A shallow, wide survey 
of 10^ SNe over 20,000 deg^ (e.g., as expected from 
LSST) could potentially limit Sxh < 3 x 10"'^ if Ux ~ 
—3, and SxH < 2 x 10~^ if Ux ^ 0. These limits 
are for a growth exponent of s = 2 in equation (jj]). 
If dark energy fluctuations are redshift-independent 




FIG. 1: Angular power spectrum {Ci^) of SNe luminosity, 
relative to the mean, as a function of redshift. The thin 
solid "lensing" line is the power spectrum due to gravi- 
tational lensing; the luminosity fluctuations generated by 
magnification associated with foreground mass perturba- 
tions. We show three examples of power spectra due to 
large-scale dark energy inhomogeneities, characterized by 
(see eq. ©) 51h = 10"'^, fcc = 0.1 h Mpc"\ s = 2 and 
Ux — —2 (bottom solid line), -I (middle dot-dashed line), 
and (top dashed fine). The two dashed "shot- noise" 
curves are for surveys with 3,000 SNe in 10 deg^ (bottom) 
and 10*^ SNe in 20,000 deg^ (top). The two stair-step fines 
show the total noise in measuring luminosity anisotropies 
(equation 6), including the contribution from lensing term 
Cf. The width of the step is equal to the bin width A£. 



(s = 0), then the limits improve by a factor of 2 to 
3. Independent of Ux and s, LSST could limit the 
dark energy quadrupole (/ = 2) to a fractional energy 
density as low as 2 x 10~^. 

The phenomenological dark energy fluctuation 
power spectrum may be observationally determined 
by mapping the luminosity residuals as a function of 
SN redshift, and then measuring the angular power 
spectrum of the residuals (i.e., in the same manner 
that CMB residuals are treated). The imprint of 
dark energy perturbations on the CMB, and cross- 
correlations with other tracers of large scale structure, 
is limited by cosmic variance. By using a shallow, all- 
sky SN survey the dark energy perturbation signal is 
dramatically improved. Thus, beyond the CMB, it 
is likely that SN data will become the most powerful 
probe of anisotropic expansion and dark energy inho- 
mogeneity in the universe. 

The multipole moments of the luminosity residu- 
als, (<5L)fm, can be used for cross-correlation stud- 
ies with other trace fields of the large-scale structure. 
The moments can be used for additional statistical 
studies to constrain signatures of global anisotropies, 
which may be due to anisotropic dark energy or other 



4 



direction-dependent phenomena. In a similar fashion, 
such studies wiU benefit from an all-sky SN sample, 
which could be obtained with a shallow SN survey 
out to z ~ 1 rather than a deep (but narrower) survey 
that attempts to recover SNe at z ~ 2. Alternatively, 
the absence of SN luminosity fluctuations due to dark 
energy inhomogeneities consists of a null test of the 
cosmological constant. 

The cosmological community is aggressively devel- 
oping and perfecting probes of the time dependence 
of the dark energy equation-of-state. It is to be noted 
that time-dependent models of the dark energy also 
necessarily predict spatial dependences. Type la su- 
pernovae have become one of the most powerful probes 
of the dark energy equation of state, and are expected 
to provide precision constraints on possible time vari- 
ation. We show that these same supernovae also offer 
one of the most powerful methods to constrain spa- 
tially inhomogeneous models of the dark energy. 

We have presented a simple phenomenological 
model for large-scale spatial fluctuations of the dark 



energy, and have derived limits to how well this model 
can be constrained with various proposed surveys. It 
will be useful to quantify the effective variations in 
dark energy for a variety of theoretical models. In 
addition, it would be worthwhile to consider other 
potential cosmological probes of dark energy inho- 
mogcneity, including signatures in weak lensing shear 
measurements of galaxy shapes, as well as spatial in- 
homogeneities in baryon acoustic oscillation measure- 
ments. The advent of almost all-sky SN samples, such 
as would be expected from LSST, will sharply con- 
strain the amount of spatial variation in the dark en- 
ergy. 
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